Dyson- Schwinger approach to color superconductivity at finite temperature and 

density 
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We investigate the phases of dense QCD matter at finite temperature with Dyson-Schwinger 
equations for the quark propagator for Nf = 2 + 1 flavors. For the gluon propagator we take a 
fit to quenched lattice data and add quark-loop effects perturbatively in a hard-thermal-loop-hard- 
dense-loop approximation. We consider 2SC and CFL-like pairing with chiral up and down quarks 
and massive strange quarks and present results for the condensates and the phase diagram. We find 
a dominant CFL phase at chemical potentials larger than 500 — 600 MeV. At lower values of the 
chemical potential we find a 2SC phase, which also exists in a small band at higher temperatures 
for larger chemical potentials. With values of 20 — 30 MeV, the critical temperatures to the normal 
phase turn out to be quite small. 



I. INTRODUCTION 

The phase diagram of strong-interaction matter (QCD 
matter) is a central object in experimental and theoret- 
ical studies [HQ- At low temperatures and densities, 
quarks and gluons are confined into hadrons while at high 
temperatures free quarks and gluons are the dominant 
degrees of freedom. 

The low-density regime of QCD matter is understood 
quite well as it is accessible by heavy-ion collisions and, 
on the theory side, by first-principles calculations on the 
lattice. These lattice QCD simulations at physical quark 
masses show a crossover transition from the hadronic 
phase at low temperatures to the quark-gluon plasma 
at high temperatures 0, H| . 

On the other hand, the regime at non-vanishing net 
baryon densities is much less understood. Lattice QCD 
fails due to the fermion sign problem that prevents Monte 
Carlo simulations to be performed at non-zero baryon 
chemical potential, and only a small window can be ac- 
cessed with extrapolation methods 0411 ■ 

The regime of large densities and low temperatures 
is totally inaccessible with lattice QCD. There, quarks 
form Cooper pairs which condense as color superconduc- 
tors PJJ-[l6| . At very high densities the coupling becomes 
small and QCD can be studied with weak-coupling meth- 
ods. In this way, it was shown that the ground state is 
a color superconductor in the color-flavor-locked (CFL) 
phase, where up, down, and strange quarks are paired 
symmetrically p7H19| - On the other hand, at quark 
chemical potentials below 1 GeV, as relevant for the de- 
scription of compact-star interiors, weak-coupling meth- 
ods are not applicable, and one has to resort to non- 
perturbative methods. In this regime the strange-quark 
mass cannot be neglected against the chemical potential, 
eventually making the CFL pairing pattern unfavorable. 
This has been studied extensively in effective models, 
like the Nambu-Jona-Lasinio (NJL) model. It was found 
that at intermediate chemical potentials the 2SC phase 
is favored [20| . in which only up and down quarks form 
Cooper pairs. Moreover, when the constraints of beta 



equilibrium and electric neutrality are taken into account 
in addition, these model studies show a rich phase struc- 
ture with many different types of color-superconducting 
condensates [2l| - t23| . 

Unfortunately, these results are very sensitive to the 
model parameters. In order to make quantitative 
predictions, it is therefore necessary to employ non- 
pcrturbative methods which are directly rooted in QCD. 
In this context, Dyson-Schwinger equations (DSEs) are 
a promising approach. They have been used to investi- 
gate the chiral and deconfinement transition a nd g ood 
agreement with lattice results was obtained [13, l25l | In 
Refs. [H, [27]], these studies were extended to finite den- 
sity. It was found that the resulting phase diagram con- 
tains a critical endpoint beyond which the chiral phase 
transition becomes first-order. 

Some time ago, the framework of DSEs was also ap- 
plied to investigate color-superconducting phases at zero 
temperature [28l - l3G | . In contrast to the NJL- model stud- 
ies, these calculations showed that the CFL phase is the 
dominant phase down to rather low chemical potentials, 
leaving practically no room for other quark phases, like 
the 2SC phase. 

In their exact formulation, DSEs are an infinite set of 
equations, which have to be truncated in practice. This 
should be done in a way that the most important features 
of QCD are included, while keeping the numerical effort 
on a reasonable level. Obviously, the truncated DSEs are 
no longer exact and typically require some model input 
for higher-order n-point functions. However, unlike effec- 
tive models, they are applicable at all energy scales and 
the approximations are systematically amendable. 

Recently, considerable progress was made in extracting 
some of the relevant input from the lattice. Specifically, 
there exists a temperature dependent parametrization of 
the quenched gluon propagator [25[ , which, together with 
a phenomenological ansatz for the quark-gluon vertex, 
has successfully been used as input for the DSE calcula- 
tions of Ref. [Is, [27|. Compared with earlier truncations, 
these gluon interactions turn out to be more attractive, 
leading to stronger non-perturbative effects. 

One goal of the present paper is therefore to perform 
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FIG. 1. Dyson- Schwinger Equation for the full quark propa- 
gator. Plain lines represent quark propagators, the curly line 
the gluon propagator. Thick dots represent dressed quanti- 
ties. 



an update of the DSE analysis of color-superconducting 
phases presented in Ref. [29|, using this new, more re- 
alistic gluonic input. We will show that this changes 
the results qualitatively, making the existence of a 2SC 
phase at low temperatures possible. The second goal of 
this work is to extend these calculations to non-zero tem- 
perature, which was not studied in Ref. (29j . 

This paper is organized as follows. In Sections [TTI and 
IIII1 we introduce the quark DSEs and the structure of 
the quark propagators in a color superconducting envi- 
ronment. Our truncation scheme is presented in Section 
IIV1 In Section [Vj we discuss our results for the color- 
superconducting condensates and for the phase diagram. 
All calculations are performed for Nf = 2 + 1 flavors. 
Thereby, we restrict ourselves to the case of a flavor- 
independent chemical potential, postponing the analysis 
of electrically neutral matter to a later analysis. Conclu- 
sion are drawn in Section [VII 



II. 



THE QUARK DYSON-SCHWINGER 
EQUATION 



The dressed quark propagator S(p) is the solution of 
the Dyson-Schwingcr equation (DSE) 



(1) 



which is depicted diagrammatically in Fig. [T] Here 
p = (pi,p) is the Euclidean 4-momentum. To describe 
the system at temperature T and chemical potential fi 
we use the Matsubara formalism, so that p^ = u) n + i\i 
where ui n = (2n + l)wT are the fermionic Matsubara fre- 
quencies. Sq(p) denotes the bare quark propagator and 
Z2 is the quark wave function renormalization constant. 
The self-energy E(p) is given by 



Z 2 V{p) = 9 2 ^TfS(q)T b l/ (p, q )Df l/ (k 



p-q) (2) 



with the dressed gluon propagator D° u (k) and the bare 
and dressed quark-gluon vertices, gTf,' and gT b (p, q), 
respectively, a and b denote color indices, \x and v 
Dirac indices and g represents the QCD coupling con- 
stant. Furthermore, we have introduced the notation 
J q = T^2 J (2tt) s > wnere the sum runs over the Mat- 
subara frequencies corresponding to (74. 



As we work in Landau gauge, the gluon propagator is 
transverse. In a thermal medium it can have two com- 
ponents in the transverse subspace, 

DfM = ^-P^M + ^-P^k), (3) 

with the 3-dimensional transverse and longitudinal pro- 
jectors 



(k) 



k^ kj 



- T^ y {k) P^ u (k) 



and the 4-dimensional transverse projector 



k^ku 



(4) 



(5) 



The gluon dressing functions Z^ b T (k) and Z^ b L (k) coin- 
cide in vacuum but, in general, differ in the medium, as 
the medium breaks the 0(4) symmetry of the vacuum. 

The quark spinors are four-dimensional objects in 
Dirac space and have 3 x Nf components in color-flavor 
space, where Nf is the number of flavors. In this arti- 
cle we generally consider Nf = 3, but sometimes refer 
to two-flavor systems for comparison. In addition we in- 
troduce the two-dimensional Nambu-Gorkov (NG) space, 
which easily allows to incorporate color superconductiv- 
ity in the formalism [3l|, HH . The NG spinors are defined 
by 



1 



V2 \^ 



V' 



1 
71 



4> T c) (6) 



with the charge conjugation matrix C = 7274- 

We take the conventions from [IJ|, [28| and parametrize 
the propagators and self-energies in NG space as follows: 



MP) 



S(P) 



_ ^0 



s+(p) 

S -(p) 



S+(p) T-(p) 
T+(p) S-( P ) 



" W{P) S"(p) 



(7) 



(8) 



(9) 



The components diagonal in NG space represent 
the normal propagators and self-energies for particles 
(+) and charge conjugate particles (— ), while the off- 
diagonal components are related to color superconduc- 
tivity. Therefore the bare propagator is diagonal in 
NG space. Without color superconductivity, also the 
dressed propagator and the self-energy are diagonal, and 
the Dyson- Schwinger system decouples into two equiv- 
alent gap equations for the quark propagator and the 
charge conjugate propagator. On the other hand, if 
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color-superconducting condensates are present, quarks 
and charge conjugate quarks are coupled. 

This becomes evident when we insert the above ex- 
pressions into the DSE, Eq. ([1}. One then obtains the 
following set of equations, which are coupled through the 
color-superconducting condensates: 



S 



±-i 



Zo s, 



(sr 1 



(10) 



By construction, the two NG-spinors Eq. ((6|) are not 
independent. As a consequence, the + and — components 
of the propagators are related to each other by 



S + (p) = -CS-(- P fC, 

r+(p) = -CT+(-p) T c, 

S + (P4,P) = 74<5 ,+ (-P4,p) t 74, 
T + (j>i,p) = 74r _ (-p 4 ,p) t 74, 



(11) 



(12) 



and analogously for the self-energies [13J, [28j 
The bare quark propagator is defined by 

(S^) _1 (p) = -ip + Z m m,f 



(13) 



with the bare quark mass m / and the mass renormaliza- 
tion constant Z m for the quark flavor /. It follows that 
(So) takes the same form with fj, replaced by — /i. 

The quark renormalization constants Zi and Z m are 
determined by the renormalization condition in vacuum 



s-\ P 2 



(—i.ft + m f )\ p 2 =u 2 



(14) 



III. PARAMETRIZATION OF 
COLOR-SUPERCONDUCTING PHASES 

In the following we restrict our calculations to spin-zero 
color-superconducting phases, which should be preferred 
over spin-one phases in most regions of the phase dia- 
gram. We are interested in the realistic case with two 
light quark flavors and a heavier strange quark. For sim- 
plicity, we consider massless up and down quarks but 
massive strange quarks. This corresponds to an interme- 
diate situation between two idealized cases: 

For infinitely heavy strange quarks, the strange sector 
decouples, and we are left with a two-flavor system. In 
this case, the 2SC phase is favored, where red and green 
up and down quarks form a diquark condensate, whereas 
blue quarks are not involved in the pairing. This leads 
to a breaking of the SU C (3) color symmetry to a SU C {2) 
subgroup, while the SUl{2)®SUr{2) chiral symmetry re- 
mains intact. The other extreme is to have three massless 
quark flavors. Here the favored phase is the CFL phase 
[1 7I | , where all colors and flavors participate in the pairing 
in a symmetric way, so that the SU C (3) color symmetry 
and the SUl{3) <£) SUr(3) chiral symmetry are broken 
down to a residual SU c +v(3) symmetry. 

Turning on a finite but small strange-quark mass, a 
CFL- like phase can still be realized. In this phase all col- 
ors and flavors arc paired in a diquark condensate, but 
there is only a residual SU c+ v(2) ® C/ c +i/(l) symmetry, 
while the full CFL symmetry is only approximate. Alter- 
natively there could be a 2SC phase, where the strange 
quarks remain unpaired. 

In order to find the solutions of the DSE in a given 
color-superconducting phase, we follow Ref. [29| and ex- 
pand the propagators and self-energies as 



s + ( P ) = Y. s t(p)Pu 

i 

T+(p) = ]£X(p)M< 



(17) 



at some arbitrary renormalization point v. 

The vertices also live in NG space. The bare vertex is 
diagonal 







r a,0 ry [ 7m 2 

V -7^ 



t =: Zif7m" 



A„ 



(15) 



where A a are the 8 Gell-Mann matrices in color space. 
Z\p is the renormalization constant of the quark-gluon 
vertex. The full vertex has off-diagonal elements in gen- 
eral, 



KM 



(p,q) A£-(p,g) 
Aj.+ (p,g) rjTfag) 



(16) 



In this work, we only take into account the diagonal func- 
tions which will be specified later, when we define our 
truncation scheme. 



and 



E+(p) = J2 E tiP)Pi, 

i 

§+(p) = 5>+(p)M„ 



(18) 



with Pi and Mj being matrices color-flavor space. These 
matrices have been constructed such that they respect 
the symmetries of the considered phase and are com- 
plete in the sense that they yield a closed set of self- 
consistency equations when inserted into Eq. (|10j) . In 
the most general case, considered here, i.e., the CFL- 
likc phase with non-vanishing strange-quark mass, both 
sets, {Pi} and {Mi}, consist of seven matrices, which are 
listed in Appendix [A] More symmetric cases, such as the 
2SC phase or non-superconducting phases, are contained 
in this parametrization as special limits, which are dis- 
cussed in Appendix [X] as well. 

The Dirac structure of the propagator is parametrized 
as 29. 33] 
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Ttip) = (74 J|T+» + l±T+.(p) + T+.(p) + |^») 75. 



(19) 



(p) accounts for chiral symmetry breaking in the nor- 
mal propagator. Df(ji) as well as Tg and T^ i (p) 
are only non-zero for color-superconducting phases with 
finite strange-quark masses. The self-energies are decom- 
posed analogously. 

Diquark condensates can be calculated as 

d = (tp T Cj 5 Oi^) = -Z 2 >/> ( l5 OiT-(q)) (20) 



where Oi is an operator in color-flavor space. In partic- 
ular, we define C uc i with 



M 2 



(21) 



describing the mutual pairing of (red and green) non- 
strange quarks in the 2SC phase as well as in the CFL 
phase, and C uc is with 



2 



H 

2 



A7 
2 



(22) 



(M 6 + M 7 - M 4 - M s ) . 



describing the pairing of non-strange quarks with strange 
quarks in the CFL phase. 

These condensates represent the dominant pairing pat- 
terns in the color-superconducting phases under consid- 
eration. In a phase with exact CFL symmetry, C uc i and 
Cuds are equal, while in a CFL-likc phase with heavier 

I 



I 

strange quarks, they can be different. In the 2SC phase, 
finally, C u d s vanishes. 



IV. TRUNCATION SCHEME 

The quark DSE (TTJ and all other equations discussed so 
far are exact QCD equations. However, in order to solve 
them, we need to specify the dressed gluon propagator 
and the dressed quark-gluon vertex, which enter the self- 
energy, Eq. In principle, they are given by their own 
DSEs but, since these involve even higher n-point func- 
tions, an exact solution would result in an infinite tower 
of equations. Therefore we have to rely on truncations. 
In this section we specify the truncation scheme used in 
the present work. 



A. Dressed quark-gluon vertex 

We begin with the dressed quark-gluon vertex, 
Eq. ([To]) . We restrict it to the diagonal NG components 
and an abclian ansatz 



i>,<?) = ^r(fc), 



(23) 



with A a as defined in Eq. (|15p and a dressing function T, 
which only depends on the gluon momentum k = p — q. 
For this function, we take a model ansatz, similar to the 
vertex proposed in Ref . [13 : 



T(k) 



Z2Z3 



di 



A 2 



Poa{y) ln(fc 2 /A 2 + l) 

47T 



(24) 



Here Z% is the ghost propagator rcnormalization con- 
stant. The second term in parentheses describes the per- 
turbative running in the ultraviolet, with /3o = (WN C — 
2N f )/3, 5 = ~9N C /(UN C - 8N f ) and a scale factor A. 
a{v) = g 2 /47r is the strong fine-structure constant at 
the renormalization scale v. Following Ref. [24[, we take 
a{y) = 0.3 and A = 1.4 GcV. For the first term, which 
is relevant for the infrared behavior, we take di = 0.5 
GeV 2 [HI and fit d\ to obtain a critical temperature 
of around 150 MeV for the chiral phase transition at 
/x = 0. We get a value of d\ = 9.6 GcV 2 . The vertex 
in Ref. [24| has additional contributions motivated by a 



I 

Slavnov- Taylor identity and a Ball-Chiu vertex construc- 
tion [34)]. As there is no strict argument for the neces- 
sity of these components, we dropped them as they lead 
to instabilities in the iteration of color-superconducting 
phases. 



In Landau gauge, the renormalization constant Z\p, 
which appears in the bare vertex Eq. (|15[) , can be ex- 
pressed as Zip = Z 2 / 'Z 3 using Slavnov- Taylor identities. 
The normal and anomalous components of the self-energy 
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FIG. 2. DSE for the gluon propagator. Curly, dotted and 
plain lines represent gluon, ghost, and quark propagators, re- 
spectively. Thick dots indicate dressed quantities. 



FIG. 3. Unquenching of the gluon propagator. The shaded 
propagator is the full Yang-Mills gluon propagator. 



B. Dressed gluon propagator 



The gluon DSE is depicted in Fig. [5] Restricting this 
equation to the first three lines corresponds to the pure 
("quenched") Yang-Mills system, while the last diagram 
describes the coupling to the quarks. Solving the Yang- 
Mills equations is numerically demanding already with- 
out quarks [35rf37j . In recent years significant progress 
has been made in this sector by combining continuum 
methods with lattice calculations. Therefore we con- 
sider a truncated version of the full gluon DSE, where 
we take the solution of the quenched Yang-Mills equa- 
tions as input and include the quark effects only per- 
turbativcly. This scheme is depicted in Fig. [3] It ne- 
glects back-coupling effects of the quark propagator on 
the Yang-Mills sector but should carry the dominant ef- 
fects of the quark loop on the gluon propagator. 

The gluon DSE is then given by 



D, 



-l.ab 



(k) = D 



-l,ab 
tiv,YM 



(k) + K b M 



(26) 



Eq. ([2]) then become 
E+(p) = 4 7 ra(,)^^7^5+(g) 7 ^<(fc), 



$+(p) 



4- z 3 

q 

4 ™(^^f r*(g)7,yJ#(*), 

(25) 



where k is again the gluon momentum p ~ q. Note that 
Z3 drops out when Eq. (|24|l is inserted. 



with the Yang- Mills gluon propagator D°^ vYM and the 
gluon polarization tensor IK* . The former has the same 
transverse Dirac structure as the full gluon propagator, 
Eq. ([3]), and is diagonal in color space, 



D^ym^) 



(27) 

For the dressing functions we adopt the paramctrizations 
of Ref. IH, 



7 YM 
J TT,TL 



(*) 



2 A 2 



k 2 A 



(fc 2 + A 2 ) : 



k 2 + qt L A 2 



bT - L ( My) ln(fc 2 /A 2 + lp 

+ A 2 \ Air 



(28) 



which have been obtained by fitting lattice results. The 
scale factor A and the parameters describing the pertur- 
bative behavior in the UV are the same as in Eq. ([24)) . 
The exponent 7 = (-13iV c + 4JV/)/(22JV c - 4N f ) is 
related to 8 by 28 + 7 = —I. The infrared behavior 
is parametrized by a temperature-independent constant 
c = 11.5 GeV 2 and temperature-dependent parameters 
qt,l and br,L, which are tabulated in Ref. [25j for var- 
ious temperatures. We determine their values at other 
temperatures by linear interpolation. 



The polarization tensor is given by 

<(fc) = -y^Tr (rfS(p)T b l ,(p, q )S( q )) (29) 
1 

with p = k + q. The trace has to be carried out in Dirac, 
color, flavor and Nambu-Gorkov space. As quantum cor- 
rections must not change the transverse nature of the 
gluon, see Eqs. ([SJ and (|2"7| . we require the polarization 
tensor to be transverse as well, i.e., 



n#(*0 = nf T (fc)PJ y (fc) + n^(fc)F^(fc). (30) 
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The dressed gluon propagator is then given by 



Z TT( k ) pT ( , \ , Z TL( k ) 

' ^ 1 Z™(k)nf T {k) » vK > k 2 + Z™{k)Uf L (k) "" W " 



Since the gluon polarization tensor depends on the 
dressed quark propagator, see Eq. (|29| . this equation 
must, in principle, be solved self-consistently together 
with the quark DSE. However, in a first step, we per- 
form a simpler non-self-consistent approximation, which 
was also employed in Ref . [IH Ull . 

C. HTL-HDL approximation 

In this scheme, we neglect the vacuum part of the po- 
larization loop and calculate the medium corrections in 



hard-thermal-loop-hard-dense-loop (HTL-HDL) approx- 
imation, using the bare propagator of massless quarks, 

(p) = -i(w„ + zm)74 - if. (32) 

In the HTL-HDL approximation, it is assumed that the 
external momenta are small compared to the chemical 
potential or temperature. This allows for an approximate 
analytical evaluation of the loop integral. The result can 
be found in textbooks, e.g., Ref. [HI, and is given by 



Z^(k)U?r b (k) = m 2 TT 
Z™( k )Uf{k) = 2< w ' 



\k\ 



k 2 



%Q 



1 —iC 



' w 

sab 



rah 



(33) 



1*1 V l fc l 



with 



D. Strange-quark mass 



iQ(ix) 



1 ix+1 
— In 

2 ix-1 



arctan 



(34) 



and bosonic Matsubara frequencies uj m = 2rmrT. The 
effective transverse and longitudinal gluon masses are de- 
fined by 



m TT,TL = NfO>TT,TL{k) 



ttT 2 



— ) 

IT ' 



with the renormalization-point independent running cou- 
plings, which are given by [28j 



OtTT,Th{k) 



T(k)Z^^ TL (k) 
Z2Z3 



a(v). 



(36) 



In this approximation, all quark-mass effects and color- 
superconducting contributions to the polarization tensor 
are neglected, but, as it provides analytical results for 
the gluon polarization, it has the advantage to keep the 
numerical effort at the same level as in a pure rainbow 
truncation with quenched gluons. 

An extended truncation scheme, where the solutions 
of the quark DSE are self-consistently used in the gluon 
polarization tensor, will be discussed in a forthcoming 
paper (39| . 



The last input quantity to be specified is the strange- 
quark mass at the renormalization scale v. Accord- 
ing to the particle data group (PDG) [4(|, its value is 
m s = 95 ± 5 MeV in the MS renormalization scheme at 
a renormalization scale of v = 2 GeV. In earlier Dyson- 
Schwinger calculations [29j | the renormalization point was 
therefore chosen to be v — 2 GeV as well, and the PDG 
value for m s was directly used as an input. 

However, as already pointed out in the Introduction, 
the gluon dressing functions and vertices used in Ref. [29| 
differ from ours, which has a strong effect on the mass 
functions. This is illustrated in Fig. 21 where the vacuum- 
mass functions of the (chiral) up and down quarks are 
shown for the two interactions used in Ref. [291 ] and for 
the present setting. 

An essential difference is that, for the former, chiral 
symmetry restoration sets in at lower momenta. As a 
consequence, p = 2 GeV already belongs to the pertur- 
bative regime, so that this scale can be used as renormal- 
ization point. 

In contrast, with the improved parametrizations of 
quark-gluon vertex and gluon propagator used in our cal- 
culations, there are still considerable non-perturbative ef- 
fects at 2 GeV, and therefore we have to renormalizc at 
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FIG. 4. Vacuum-mass functions M{p) = B+{p)/AUp) of 
chiral quarks (m = 0) for the interactions oliji from [2t| and 
in our setting. 



a higher scale to be in the perturbative region. To be 
on the safe side, we choose v = 100 GcV. Thus, in order 
to make contact to the PDG value of the strange-quark 
mass, we have to evolve it to that higher scale, employ- 
ing the same perturbative running which was used in the 
PDG analysis. 

The scale evolution of the quark mass is then described 
by the differential equation 



, dm(y) 



while the coupling obeys the equation 
2 da(v) 



dv 2 



(37) 



(38) 



where the functions /3(a) and 7(a) are known to four- 
loop order for MS renormalization [4(| SJ. With the 
initial conditions m s (2 GeV) = 95 MeV and a(Mz = 
91.2 GcV) = 0.118 we integrate Eq. (37]) to m b = 4.2 
GeV with 4 flavors and finally to 100 GeV with 5 flavors, 
giving m^ fS (100 GeV) = 54 GcV. This strange-quark 
mass was also chosen in [27j , although at a slightly lower 
renormalization scale of 80 GeV. 

So far, we have neglected the fact that the present 
renormalization scheme corresponds to a momentum sub- 
traction (MOM) scheme, which is different from MS. To 
estimate the effect, we consider the leading-order running 
coupling, which is given by 



a{v) 



4tt 



{ll-2/W f )\og{p*/k% CD ) 



(39) 



with the QCD scale Aqcd = 213 MeV for 5 flavors and 
296 MeV for 4 flavors in MS renormalization [H . This 
leads to a one-loop mass scaling 



m(y) = m(va) 



l°g^ 2 /A| CZ5 ) 



(40) 



FIG. 5. Vacuum-mass functions M(p) of chiral quarks (m = 
0) and strange quarks (m = 30 MeV, m — 54 MeV). 



with the anomalous mass dimension 7 = 12/ (33 — 27V/). 
The MOM scheme has a bigger scale, which can be re- 
lated to the M£-scale as A-^ = 0.346A M om 0, S3- 
Using this scale to evolve the quark mass to 100 GeV 
with Eq. PO"]) and rescaling the result with the ratio be- 
tween the 4th-order and leading-order result for MS scal- 
ing to roughly include higher-order corrections, we obtain 
a mass of around m s (100 GcV) = 35 MeV. 

We therefore expect the masses in the MOM scheme 
to be lower than in MS and perform calculations for 
m s = 54 MeV and m s = 30 MeV at our renormalization 
point of 100 GeV. These values can be seen as upper and 
lower bounds for the strange-quark mass. 

The resulting vacuum-mass functions, together with 
the corresponding function for chiral quarks are shown 
in Fig. [5] The mass functions seem to be unrealistically 
large, when we compare the p = values with typical con- 
stituent quark masses in phenomenological models. In an 
improved truncation scheme, which will be discussed in a 
future publication [39j . we obtain smaller masses, which 
are of the order of 450 MeV for up and down quarks and 
600 to 700 MeV for strange quarks. It should be empha- 
sized, however, that the quark masses are not observable. 
Moreover, the size of the non-perturbative regime will 
roughly stay the same. 



V. RESULTS 

In Fig. |S] we show the diquark condensates defined in 
Eqs. (201) " (H2) at l° w temperatures, T = 10 MeV, as 
functions of the chemical potential. As discussed in Sec- 
tion [Till our formalism allows for the description of both 
2SC and CFL-like pairing. Since the condensate C u d, re- 
lated to the mutual pairing of up and down quarks, is 
present in both phases, we indicate explicitly in which 
phase the solutions have been obtained. We also do 
so for the condensate C u d s , although it is non-zero only 
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FIG. 6. Dependence of 2SC and CFL condensates at T = 
10 MeV on the chemical potential for m s = 30 MeV (top) 
and m s = 54 MeV (bottom). All condensates are given in 
arbitrary, but equal, units. 



in the CFL phase. In addition, there is always a non- 
superconducting solution, where all diquark condensates 
vanish. This is not shown in the figure. 

The upper and the lower panel of Fig. [5] correspond to 
m s = 30 MeV and 54 MeV, respectively. The qualitative 
behavior is similar in both cases, and even the quantita- 
tive differences are not large. All condensates rise with 
increasing chemical potential. At lower values of /i, we 
only find a 2SC solution of the DSE, whereas above a 
threshold of about 500 or 600 MeV, depending on m s , 
there is also a CFL-likc solution. 

This behavior is related to the //-dependence of the 
dynamical quark masses, which are displayed in Fig. [7] 
Since the normal and color-superconducting dressing 
functions are mutually coupled through the DSEs, the 
dynamical quark masses depend on the phase as well. 
At low \i we find a non-superconducting solution with 
spontaneously broken chiral symmetry, where the up and 
down quarks have relatively large masses (dash-dotted) 
line. On the other hand, there is always a 2SC solution 
with massless up and down quarks. In the CFL-like so- 



FIG. 7. Dynamical quark masses M(0) = B+(0)/C + (0) 
as functions of the chemical potential at T — 10 MeV for 
m s = 30 MeV (top) and m s = 54 MeV (bottom). The 
up- and down-quark masses are only shown for the non- 
superconducting chirally broken phase. They vanish in the 
2SC phase and are very small in the CFL phase. 



lutions, the up and down quarks are almost massless as 
welQ. The strange quarks, on the other hand, are quite 
heavy at low and intermediate chemical potentials. As 
a consequence, the Fermi momenta of strange and non- 
strange quarks arc very different, so that their mutual 
pairing is inhibited. However, with increasing chemical 
potential, the strange quarks become lighter and even- 
tually undergo a crossover to even smaller masses. This 
effect, which is seen in the 2SC phase (dashed line), also 



1 In the CFL phase, the up and down quarks have a small mass 
(< 2 MeV), which is induced by their coupling to the massive 
strange quarks. This is possible because the CFL condensates 
break the chiral SU (2) symmetry. The situation is further com- 
plicated by the fact that red up, green down, and blue strange 
quarks are mixed in the CFL phase (see Eq. dAll ll, In princi- 
ple, one should therefore project on the different quasi-particlc 
modes. However, for simplicity, we just show the Pg component 
in Fig. [7] corresponding to red and green strange quarks. 
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FIG. 8. Temperature dependence of 2SC and CFL conden- 
sates at n = 580 MeV for m s = 30 MeV (top) and at /u = 680 
MeV for m s = 54 MeV (bottom). The units are the same as 
in Fig. E] 



triggers the onset of CFL-like solutions. In the latter the 
strange-quark mass (dotted line) is slightly lower than in 
the former, but the difference is small and almost van- 
ishes at very high /i, where the mass is dominated by 
pcrturbativc effects. 

As discussed earlier, the quark masses are neglected in 
the gluon polarization function in the HTL-HDL trunca- 
tion. Therefore, if there is no pairing between strange and 
non-strange quarks, the strange sector decouples from the 
non-strange sector. As a consequence, the strange-quark 
mass is not affected by a possible chiral phase transition 
of the up and down quarks. Likewise, C^f c does not de- 
pend on m s and is not influenced by the rapid change of 
M s in the crossover region. 

In Fig. [8] the temperature dependence of the diquark 
condensates is shown for m s = 30 MeV at a chemi- 
cal potential of /a = 580 MeV and for m s = 54 MeV 
at fi = 680 MeV (lower panel). The 2SC conden- 
sate smoothly decreases with increasing temperature and 
eventually vanishes at a temperature between 20 and 
25 MeV. In contrast, the CFL-like solution ceases to 



exist already above a lower critical temperature, where 
the condensates are still finite. This behavior suggests 
that, with increasing T, the system first undergoes a first- 
order phase transition from the CFL phase to the 2SC 
phase, followed by a second-order phase transition to the 
normal-conducting phase. 

However, if there are several solutions of the DSE at 
given temperature and chemical potential, the favored 
phase can only be determined by comparing the corre- 
sponding pressures. The pressure is equal to the effective 
action, which in general is given by 

r = ^TrlnS-^-^Tr (l - Z 2 S^(p)S(p))+T 2 . 

p p 

(41) 

Here T 2 denotes the two-particle irreducible part. Our 
present truncation with HTL-HDL approximation corre- 
sponds to the expression 
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(42) 

The resulting pressure is quartically divergent. Although 
finite expressions can be obtained by calculating pressure 
differences between competing phases [29j , unfortunately 
the results are numerically very unstable. This is due 
to the fact that the convergence behavior is determined 
by the high-energy tails of the integrands, which contain 
only little information about the phase under considera- 
tion. 

We therefore proceed as follows. As argued in ap- 
pendix |Bj there are reasons to believe that the solu- 
tions found by solving the DSEs iteratively, correspond 
to global or local maxima of the pressure, i.e., to stable or 
metastable solutions, but not to unstable ones. Accord- 
ing to this hypothesis (which is often assumed in DSE 
calculations), this means that we always find the correct 
solution in the case of second-order phase transitions. 

For first-order phase transitions, the situation is more 
complicated. In this case, there is a certain regime, where 
a stable and a metastable solution of the DSE exist, 
which both can be found by iteration. The exact po- 
sition of the phase transition, which manifests itself in 
a jump between the two solutions, can then only be de- 
termined by studying the pressure of the system. How- 
ever, because of the numerical uncertainties mentioned 
above, in practice this does not further narrow down the 
phase-transition region. We therefore restrict ourselves 
to calculating the spinodal lines of the first-order region, 
i.e., the lines where the metastable solutions disappear. 
As these regions mostly have only a small extent, they 
still give a good estimate for the phase boundary. 

The resulting phase diagrams are shown in Fig. [§] 
Again, the upper and lower panels correspond to m s = 30 
MeV and 54 MeV, respectively. As pointed out ear- 
lier, in the HTL-HDL approximation, the strange and 
non-strange sectors decouple, except for the CFL phase, 
where they arc coupled by the condensate. Therefore the 
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FIG. 9. Phase diagram for m s = 30 MeV (top) and m s — 
54 MeV (bottom). The shaded areas indicate the region of 
the first-order phase transition, bounded by spinodal lines 
(dashed). Solid lines indicate second-order phase transitions, 
CP denotes the tricritical point. 

only effect of m s on the phase diagram is a shift of the 
CFL phase boundary. 

The shaded areas indicate first-order regions bounded 
by spinodal lines. The most prominent ones correspond 
to the chiral phase transition from the chirally broken 
phase to the 2SC phase or, at somewhat higher tempera- 
tures, to the non-superconducting chirally restored phase 
(with respect to the up and down quarks). The spinodal 
region and, hence, the first-order phase transition inside 
this region end at a tricritical point at T = 60 MeV and 
fj, = 300 MeV. Above this point, the chiral phase transi- 
tion is of second order. 

The temperature of the tricritical point is lower than 
in the DSE analysis of Ref . |2f| , where it was found at 
T = 95 MeV and n = 280 MeV. The authors use a similar 
truncation scheme but consider only two flavors. This 
also leads to a higher critical temperature of 180 MeV at 
/x = 0. As shown in Ref. [23|, the inclusion of strange 
quarks reduces the temperature of the critical point and 
shifts it to a slightly higher chemical potential. This is 



consistent with our result. However, since in Ref. [27| 
a truncation without the HTL-HDL approximation was 
used, a direct comparison of the 2 + 1-flavor results with 
ours is not possible. 

The low-^( part of the phase diagram and the region 
around the critical point are mainly shown for complete- 
ness, while our focus is on the color-superconducting 
phases at low temperatures and higher chemical poten- 
tials. Here we find a 2SC phase, followed by a CFL phase. 
The latter is favored above ~ 500 MeV for m s = 30 MeV 
and - 600 MeV for m s = 54 MeV. The 2SC phase also 
extends to somewhat higher temperature than the CFL 
phase (cf. Fig. [8]) and, thus, separates the latter from the 
normal-conducting chirally restored phase. The phase 
transition between 2SC and restored phase is second or- 
der and takes place around 20—30 MeV. Remarkably, the 
critical temperature is slowly decreasing with increasing 
temperature, despite the fact that the zero-temperature 
condensate increases (cf. Fig. |6|). 

As is clear from the discussion of Fig. [51 the phase tran- 
sition between 2SC and CFL phase at low temperature 
is first order. However, the spinodal region turns out to 
be very narrow, so that it is hardly visible in the phase 
diagra nfl Moreover, we have numerical indications that 
at large /z the phase transition is second order, i.e., the 
spinodal region ends again in a critical point. Such a 
behavior was also found in the NJL-model analysis of 
Ref. [2(j. For m s = 30 MeV, the critical point seems to 
be somewhere around fi = 900 MeV, but it is difficult to 
localize its position exactly. 

VI. CONCLUSIONS AND OUTLOOK 

We have investigated color-superconducting phases in 
QCD with Dyson-Schwinger equations in Landau gauge. 
For the gluon propagator we used a recent parametriza- 
tion of quenched lattice data and included quark-loop ef- 
fects perturbativcly in a hard-thermal-loop-hard-dcnsc- 
loop approximation. For the vertex we used a phe- 
nomenological ansatz constrained by the known asymp- 
totic behavior in the ultraviolet, while the infrared 
strength is fitted to obtain a chiral phase-transition tem- 
perature of 150 MeV at /j, = 0. 

We have studied 2SC and CFL-like phases for Nf = 
2 + 1 flavors with chiral up and down quarks, and massive 
strange quarks. We find a dominant CFL phase at high 
chemical potentials and low temperatures while there is 
always a thin 2SC band at higher temperature, separat- 
ing the CFL phase from the non-superconducting chirally 



2 This means that the 2SC solution is unstable almost everywhere 
where a CFL solution exists. The 2SC condensates shown in 
Figs. \E\ and [8] have been obtained by forcing the strange conden- 
sates to be exactly zero during the iteration. Otherwise, starting 
with a small non- vanishing strange condensate, the iteration only 
converges to the 2SC solution if we are outside the CFL regime 
or in the very small spinodal region. 
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restored phase. The transition to the latter takes place 
at critical temperatures around T = 20 — 30 MeV. 

Our most important result is that we also find a sta- 
ble 2SC phase at small temperatures and intermediate 
chemical potential, /x < 500 MeV. This is in strong con- 
trast to the results of Ref. [29j where the CFL phase 
was found to be favored practically all the way down 
to the hadronic phase at zero temperature. Although 
our calculations have been performed for T > 10 MeV 
only, a naive extrapolation of the phase boundaries to 
T = does not change our results qualitatively. The 
difference should therefore mainly be attributed to the 
improved gluon propagators and vertices, which we have 
used. This interaction is effectively stronger than those 
used in Ref. [2{|. As a consequence, the strange quarks 
stay heavy up to higher chemical potentials, inhibiting 
their pairing with the non-strange quarks. 

Of course, our analysis leaves room for various im- 
provements in the future. Perhaps the most unsatis- 
factory feature of the present truncation scheme is the 
HTL-HDL approximation, which neglects quark-mass 
and pairing-gap effects on the quark-polarization loop 
in the gluon propagator. Therefore, we arc currently 
working on an improved scheme, where the dressed 
Nambu-Gorkov propagators from the quark DSE are self- 
consistcntly used to evaluate the polarization loop. Pre- 
liminary results indicate that this leads to quantitative 
changes, like a reduction of the dressed quark masses, 
but does not alter the phase diagram qualitatively [39j . 

Further possible improvements concern the quark- 
gluon vertex, which is only constrained in the UV so 
far, while for the IR behavior we used a phenomcnolog- 
ical model ansatz. It is desirable to include non-abelian 
vertex structures, as found in Ref. [H|, and finally to 
implement a truncated vertex DSE self-consistently. 

In the present paper, we have considered a single, fla- 
vor independent chemical potential. These studies should 
be extended to electrically neutral matter in beta equi- 
librium, as relevant for compact stars. Such an analy- 
sis has been performed in Ref. [30j for zero temperature 
using the "old" gluonic interaction. Since electric neu- 
trality splits the Fermi surfaces of up and down quarks, 
it disfavors the 2SC phase [Hj]. Therefore, it would be 
interesting to sec, whether this phase survives when the 
improved interaction is used. 

Another interesting question is, whether the QCD 



phase diagram contains inhomogeneous phases. Phases 
with spatially varying chiral condensates have been found 
in effective models |4(| 1471 and in large- N c QCD with 
weak-coupling methods [48l |49| . It would be exciting 
to study these phases with DSEs as well. This also 
includes inhomogeneous color superconductors (50l - [5^ |. 
which may exist when electric neutrality is enforced. 

Finally, it is an important, but difficult, task to calcu- 
late the pressure in the various phases, which would allow 
for a more precise determination of the first-order phase 
boundaries. Although an analytical expression for the 
pressure is given by the effective potential, its numerical 
evaluation was so far hindered by the poor convergence of 
the integrals. A solution of this problem would constitute 
a major breakthrough. 
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Appendix A: Color-flavor structure 

The matrices Pi and Mi, which parametrize the color- 
flavor structure of the quark propagator and self-energies, 
Eqs. (jTTJ) and (fl"8|) . are constructed such that they yield a 
closed set of self-consistency equations when inserted into 
Eq. (fT0"|) . Additional restrictions come from the require- 
ment that the resulting propagator must be consistent 
with the symmetries of the considered phase. 

In the most general case considered here, the CFL-like 
pairing pattern with finite strange-quark mass, the 
residual color-flavor symmetry is SU c+ v{2) ® C/ c +y(l), 
generated by t q — with a = 1,2,3, and 8. One 
then finds that the sets {Pi} and {Mi} consist of 
seven matrices each [29[. In the color-flavor basis 
{(r, u), (g, d), (b, s), (r, d), (g, u), (r, s), (b, u), (g, s), (b, d)} 
these are given by 



(5i,r + $i,2 k,2 <5i,4 

Si,5 Si5 Si 3 



Pi 



Si,i 



(Al) 



>i,6 



$i,7 



3»,6 



12 



f5 ty i + 5 l . 2 S ii2 S lA 



Mr 



$i,7 



(A2) 



Si., 



Si,7 



where we have slightly modified the notation of Ref. [29[ 

More symmetric phases arc contained in this 
paramctrization as special cases: 

In the non-superconducting phase all anomalous self- 
energy components vanish, = 0, and the normal self- 
energy is diagonal in color and flavor, i.e., E+ — for 
i = 2,4,5. The remaining self-energy components are 
equal for up and down quarks, but can be different for 
strange quarks, i.e., E^~ = E^" and E^ = Eg". 

In the 2SC phase, only red and green up and down 
quarks are paired, described by (p^ = — </>J, whereas all 
other components of the anomalous self-energy vanish. 
Accordingly, the normal self-energy depends on whether 
it corresponds to paired (red or green) up or down quarks 
(E^), unpaired (blue) up or down quarks (E^), or to 
strange quarks (Ejj~ — Eg^), while the non-diagonal com- 
ponents vanish again. 

In the CFL phase with three massless flavors, the resid- 
ual symmetry is SU c +v(3), generated by r a — with 
a = 1, . . . , 8, and the quasi-particle spectrum consists of 
an octet and a singlet [l7| . Therefore there are only two 
independent color-flavor components of the anomalous 
self-energy [53| . ^ = 0g = $ , 4>\ — 4>\ — and 
03 = (f)^ + $2 ■ Analogous relations hold for the normal 
self-energies. 



Appendix B: A remark on the iterative procedure 

Throughout this work we use a fixed point iteration to 
obtain the self-consistent solutions of the DSE. Schemat- 
ically, this can be formulated as an iteration 



x-i+i = <f(xi) 



(Bl) 



with the fixed points x* = tp(x*), which correspond to 
the solutions of the DSE. 

In general, however, not every solution can be found 
in this way, but the iteration only converges if it is con- 
tracting, i.e., if 



\Xi+l 



\ip{x z ) - ip{x*)\ ^ L 



\Xi 



(B2) 



with a constant L < 1. This expression is just the dis- 
cretized derivative |y>'(x*)|. 

The phase diagrams and in particular the spinodal 
lines presented in Sec. [V] are based on the assumption 
that the numerically stable solutions of the DSE, i.e., 
those solutions for which the iteration converges, also 
correspond to thcrmodynamically stable or metastablc 
solutions, i.e., to global or local minima of the thermo- 
dynamic potential. 

This conjecture is motivated by the observation that it 
holds in the NJL model, where it can be shown analyti- 
cally. To this end, we consider the standard NJL model, 
defined by the Lagrangian [|54| 

C NJL = 4>{i$ - m)il> + G [(44>) 2 + (i'*l 5 Ti') 2 ] (B3) 

with bare quark mass m and a positive four-point cou- 
pling constant G. The NJL mean-field thermodynamic 
potential in vacuum is given by (see, e.g., Ref. [14| ) 



n 



(M - mf 
AG 



2N c N f 



d 3 q 
(2ny 



■^q 2 + M 2 , (B4) 



which depends on the dressed quark mass M. Addition- 
ally a regularization needs to be specified, as the integral 
is divergent. 

The actual value of M is determined by minimizing il. 
This leads to the stationarity condition 



sn 

SM 



M 



<P(M) j. Q 



2G 



with 



<p(M) =m + AGN c N f 



d 3 c 



M 



(27T) 3 J? 



M 2 



(B5) 



(B6) 



Obviously, Eq. (|B5P is equivalent to the DSE ( "gap equa- 
tion") M = ip(M), which has the solutions M* . The 
physical stability of these solutions can then be checked 
by investigating the second derivative 



s 2 n 



SM 2 



M=M* 



1 

2G 1 



W(M*)) 



(B7) 



3 Our Pa and Afg correspond to P7 and M7 in Ref. |29| and our P7 
and Mi correspond to Pg and Mg, while there are no projectors 
Pe, and Me in that reference. 



with ip'(M) = dip J^P ■ This means, the solution corre- 
sponds to a minimum if ip'(M*) < 1 and to a maximum 
if <p'(M*) > 1. 
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Taking the derivative of Eq. (|B6[) . one obtains 



<p'(M) = 4GN c N f 



d 3 q 



(B8) 



(2tt) 3 (q 2 + A/ 2 ) 3 /2 
which is positive for all values of M. Therefore, we have 
< 1 for a minimum of f2 



W{M*)\ 



> 1 for a maximum of 



(B9) 



Comparing this with Eq. ()B2|) . we see immediately that 
the maxima of the potential correspond to numerical un- 



stable iterative solutions, while minima are numerically 
stable. This allows in principle to find all minima of the 
potential by iterating the gap equation. It can also be 
checked easily that the introduction of finite temperature 
and chemical potential does not change these results. 

For QCD DSEs, the situation is much more compli- 
cated as we have functional derivatives and an in prin- 
ciple infinite dimensional system. Therefore we are not 
able to provide an analogous analytic argument for the 
relation between physical and numerical stability. How- 
ever, we expect this relation still to be valid. 
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